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THEORETICAL  MODELLING  OF  CELL  AND  BATTERY  CHARACTERISTICS  AT  SHORT 

TIMES/HIGH  FREQUENCIES 


M.L.  Daroux,  B.D.  Cahan,  E.B.  Yeager 


Case  Center  for  Electrochemical  Sciences  and  Chemistry  Department, 
Case  Western  Reserve  University, 

Cleveland,  Ohio,  44106 


INTRODUCTION  : 

The  design  of  batteries  intended  to  deliver  very  high  current 
pulses  with  rise  times  on  the  order  of  microseconds  requires 
consideration  of  factors  usually  neglected  in  conventional  battery 
design.  In  previous  work  [1]  we  have  shown  that  the  short  time/high 
frequency  response  of  a  physically  large  battery  may  be  significantly 
limited  by  geometric  and  physical  factors.  In  constructing  theoretical 
models  which  can  be  used  to  evaluate  design  factors  in  pulse  batteries 
it  is  therefore  necessary  to  take  into  account  the  whole  battery, 
including  the  plate  geometry,  the  internal  and  external  connections, 
and  eventually  the  load  itself.  On  the  basis  of  a  modelling  study  which 
takes  some  of  these  factors  into  consideration,  this  paper  describes 
the  design  implications  of  geometric  and  physical  variables  such  as 
electrode  and  electrolyte  resistivity,  electrolyte  thickness,  and  cell 
length . 

The  details  of  the  modelling  procedure  are  described 
elsewhere  [1,2],  and  are  only  summarised  here.  The  results  described  in 
this  paper  have  been  obtained  for  a  single  cell  of  semi- infinite 
rectangular  parallel-plate  geometry  (a  stripline  cell)  with  infinitely 
thick  solid  electrodes  (Figure  1).  The  assumption  of  infinite  electrode 
thickness  does  not  affect  the  results  at  the  high  frequencies  of 
greatest  interest  since  nearly  all  the  current  is  carried  at  the 
surface  of  the  conductor,  but  will  lead  to  errors  at  low  frequencies. 
The  modelling  has  been  carried  out  by  considering  the  cell  behavior  as 
that  of  a  distributed  network  comprised  of  differential  elements  having 
the  equivalent  circuit  shown  also  in  Figure  1.  The  series  components 
are  the  electrode  inductance,  Ls,  and  the  skin  effect  resistance,  r^. 
The  parallel  components  comprise  three  series  combinations  of  a 
resistor  in  parallel  with  a  capacitor.  These  represent  the  positive 
electrode  interface  (subscript  W) ,  the  electrolyte  (subscript  E) ,  and 
the  negative  electrode  interface,  (subscript  C) ,  respectively.  The 
values  of  Ls,  rgfc  and  the  electrolyte  components  can  be  estimated  from 
the  physical  and  geometric  properties  of  the  cell  as  described 
elsewhere  [1].  For  the  calculations  reported  here  the  value  of  the 


Faradaic  resistance  of  the  negative  electrode  has  been  assumed 
infinite,  so  that  this  electrode  behaves  as  a  bare  metal  showing 
capacitative  behavior  only.  On  the  basis  of  this  equivalent  circuit  the 
impedance-  of  the  cell  has  been  calculated  over  a  wide  range  of 
frequencies  from  10^  to  10**  Hz.  By  use  of  inverse  Fourier 
transformations  it  has  also  been  possible  to  calculate  the  transient 
response  of  the  cell,  for  example  the  transient  voltage  response  to  a 
current  pulse  of  unit  amplitude. 


RESULTS  AND  DISCUSSION  : 

The  effects  of  a  number  of  the  design  variables  on  the  cell 
characteristics  have  been  calculated  using  the  approach  described 
above.  The  variables  studied  include  electrode  and  electrolyte 
resistivity,  electrolyte  thickness,  interfacial  resistances  and 
capacitances,  and  cell  length. 

Fig  2,  which  has  been  taken  from  [1],  illustrates  a  number  of 
the  different  types  of  limiting  behavior  that  can  be  calculated.  The 
results  are  presented  as  Bode  plots  for  log  Z  (the  modulus  of  the 
complex  impedance) , and  the  phase  angle  theta.  A  set  of  plots  is  given 
for  electrolyte  resistivities  ranging  from  10' 2  to  10-*  O.cm.  Over 
limited  frequency  ranges,  the  behavior  can  be  described  by  simplified 
versions  of  the  equivalent  circuit  given  in  Fig.  1.  The  illustrations 
of  this  given  at  the  top  of  Fig.  2  correspond  to  the  case  where  the 
electrolyte  resistivity  is  10^  O.cm.  For  this  particular  curve,  it  can 
be  seen  that  at  high  frequencies  (region  E)  the  cell  impedance  is 
dominated  by  the  series  inductance  of  the  electrodes  and  the  parallel 
capacitance  of  the  electrolyte.  (Note  that  no  interfacial  components 
contribute) .  LC  network  behavior  is  seen  and  the  impedance  has  the 
characteristics  of  a  pure  resistor  (log  Z  constant  and  a  phase  angle  of 
zero).  As  the  frequency  decreases  (region  D) ,  rg  becomes  less  than 
1/wCg  and  the  behavior  becomes  that  of  an  LR  circuit.  As  the  frequency 
continues  to  decrease  wLs  becomes  comparable  to  rsk  (region  C) ,  and 
once  rsk  becomes  the  dominant  series  component  the  behavior  becomes 
purely  resistive.  At  the  lowest  frequencies  shown  (region  A)  the 
interfacial  capacitance  becomes  the  dominant  parallel  term  and  the 
behavior  becomes  that  of  an  RC  network. 

Also  shown  in  Fig.  2  is  the  dependence  of  the  penetration 
length  1/alpha  on  frequency.  For  a  given  frequency,  all  but  1/e  of  the 
total  current  is  drawn  from  a  distance  of  1/alpha  from  the  terminals 
[1].  The  penetration  length  therefore  gives  a  measure  of  the  fraction 
of  the  cell  that  can  be  utilised  at  a  particular  frequency.  As  Fig.  2 
shows,  penetration  decreases  at  higher  frequencies/shorter  times.  Over 
the  frequency  range  shown,  increasing  the  electrolyte  resistivity  leads 
to  an  increase  in  the  penetration  length. 

Figure  3  shows  the  effect  of  the  Faradaic  impedance  of  the 
positive  electrode,  rp,  on  the  impedance  of  an  infinite  line.  Low 
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values  of  rp  have  the  greatest  effect,  acting  to  short  out  Cy  and 
thereby  reduce  the  overall  impedance  at  intermediate  frequencies  to  the 
value  determined  by  the  other  parallel  components.  At  very  high 
frequencies  rg  and  Ls  dominate  and  RL  line  behavior  is  seen,  while  at 
low  frequencies  rg  and  Cq  dominate  and  the  behavior  is  that  of  a  RC 
line.  The  frequency  at  which  the  effect  of  ry  shorting  out  Cy  appears 
increases  as  the  value  of  rg  decreases. 

Figure  4  shows  the  results  of  the  same  calculations  but  in 
the  form  of  the  voltage  transient  response  to  a  current  step  of  unit 
amplitude.  The  method  by  which  the  conversion  from  frequency  to  time 
domain  is  carried  out  is  described  elsewhere  [2].  At  short  times  the 
cell  behavior  is  primarily  inductive  and  the  voltage  decreases  rapidly. 
For  very  low  values  of  Faradaic  resistance  the  cell  behavior  then 
becomes  resistive  and  the  voltage  becomes  constant  with  time.  Beyond  a 
few  hundred  microseconds  in  the  low  rp  case  the  cell  becomes  inductive 
again  (voltage  decreasing  with  time)  and  this  can  also  be  seen, 
although  less  clearly,  in  the  corresponding  curves  in  the  frequency 
domain  in  Fig.  3.  At  longer  times  still,  though,  the  cell  becomes 
resistive,  and  then  finally  shows  capacitative  behavior  as  Cy  becomes 
dominant.  Similar  progressions  in  characteristic  behavior  with  time  can 
be  seen  with  the  higher  values  of  rp,  but  displaced  to  longer  times. 
Thus  for  rp  greater  than  10*  ^  Ohms  the  behavior  is  still  inductive 
after  a  millisecond. 

The  previous  examples  have  treated  cases  where  the  cell  is 
infinitely  long.  This  corresponds  in  practical  terms  to  situations 
where  the  physical  dimensions  of  the  cell  exceed  about  three  times  the 
penetration  length  a£  the  frequencies  considered.  Figure  5  illustrates 
in  Bode  plot  form  the  case  where  the  cell  is  of  finite  length,  that  is 
where  the  pulse  has  time  to  penetrate  to  the  end  of  the  cell.  At  high 
frequencies  (corresponding  to  periods  too  short  for  the  pulse  to  reach 
the  end  of  the  cell)  the  behavior  is  identical  to  that  of  the 
infinitely  long  line.  At  low  frequencies  the  cell  behavior  approaches 
that  of  a  pure  capacitor;  the  impedance  becomes  inversely  proportional 
to  frequency.  This  corresponds  to  the  discharge  of  the  cell  in  the 
conventional  manner.  In  the  intermediate  frequency  region  where  the 
length  of  the  cell  is  of  the  same  order  as  the  penetration  length, 
oscillatory  behavior  in  both  the  magnitude  of  the  impedance  and  the 
phase  angle  is  seen,  and  this  is  believed  to  result  from  constructive 
and  destructive  interference  as  the  pulse  is  reflected  from  the  end  of 
the  cell.  The  frequency  at  which  the  finite  length  of  the  cell  affects 
the  output  increases  as  expected  as  the  length  of  the  cell  decreases. 
For  practical  purposes  it  is  desirable  to  construct  a  cell  that  is  just 
long  enough  to  behave  as  if  it  were  infinitely  long  at  the  output 
frequency/pulse  length  required.  Any  extra  capacity  does  not  contribute 
to  the  performance  of  the  cell  and  therefore  reduces  the  overall 
efficiency. 

The  modelling  approach  illustrated  here  permits  quantitative 
calculation  of  the  performance  of  a  pulse  battery  at  high 


frequencies/short  times  where  many  intuitive  concepts  derived  from  an 
understanding  of  steady-state  behavior  do  not  apply.  It  is  not  possible 
to  list  here  all  the  calculated  effects  of  the  design  variables  that 
have  been  studied,  but  Table  1  illustrates,  for  a  particular  set  of 
conditions,  the  way  in  which  design  factors  can  influence  the 
characteristics  of  the  cell.  Many  of  these  effects,  particularly  on  the 
penetration  length,  would  be  difficult  to  predict  intutitively,  and  a 
factor  that  would  have  a  significant  effect  on  performance  at  one 
frequency  may  have  no  effect  at  another. 


SUMMARY  : 

The  effect  of  a  number  of  design  variables  on  the  performance 
of  cells  at  short  times/high  frequencies  has  been  calculated  by 
modelling  the  cell  as  a  distributed  network.  For  a  stripline  geometry, 
it  has  been  possible  to  calculate  that  at  sufficiently  short  discharge 
times  the  performance  of  the  cell  is  limited  by  physical  and  geometric 
effects,  independent  of  the  electrochemistry.  It  has  also  been  shown 
that  for  a  desired  output  frequency /pulse  length,  there  is  an  optimum 
size  of  cell,  and  that  scaling  up  beyond  this  size  produces  no 
improvement  in  performance. 
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